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Abstract 

Motivated by applications in network epidemiology, we consider the problem of deter¬ 
mining whether it is possible to delete at most k edges from a given input graph (of small 
treewidth) so that the resulting graph avoids certain subgraphs; of particular interest is the 
problem of determining whether it is possible to delete at most k edges so that the resulting 
graph has no connected component of more than h vertices, as this bounds the worst-case 
size of an epidemic. While even this special case of the problem is NP-complete in general 
(even when /i = 3), we provide evidence that many of the real-world networks of interest 
are likely to have small treewidth, and we describe an algorithm which solves the general 
problem in time '>n on an input graph having n vertices and whose treewidth is 

bounded by a fixed constant w, if each of the subgraphs we wish to avoid has at most r 
vertices. For the special case in which we wish only to ensure that no component has more 
than h vertices, we improve on this to give an algorithm running in time 0{{wh)‘^'^n), which 
we have implemented and tested on real datasets based on cattle movements. 


1 Introduction 

Network epidemiology seeks to understand the dynamics of disease spreading over a network 
or graph, and is an increasingly popular method of modelling real-world disease. The rise 
of network epidemiology corresponds to a rapid increase in the availability of contact network 
datasets that can be encoded as networks or graphs: typically, the vertices of the graph represent 
agents that can be infected and infectious, such as individual humans or animals, or appropriate 
groupings of these, such as cities, households, or farms. The edges are then the potentially 
infectious contacts between those agents. Considering the contacts within a population as the 
edges of a graph can give a large improvement in disease modelling accuracy over mass action 
models, which assume that a population is homogeneously mixing. For example, if we consider 
a sexual contact network in which the vertices are people and the edges are sexual contacts, 
the heterogeneity in contacts is very important for explaining the pattern and magnitude of an 
AIDS epidemic [T]. 

Our work has been especially motivated by the idea of controlling diseases of livestock 
by preventing disease spread over livestock trading networks. As required by European law, 
individual cattle movements between agricultural holdings in Great Britain are recorded by 
the British Cattle Movement Service (BCMS) [23]; in early 2014, this dataset contained just 
under 300 million trades and just over 133,000 agricultural holdings. For modelling disease 
spread across the British cattle industry, it is common to create vertices from farms, and edges 
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from trades of cattle between those farms: a disease incursion starting at a single farm could 
spread across this graph through animal trades, as is thought to have happened during the 
economically-damaging 2001 British foot-and-mouth disease crisis |18j . 

We are interested in controlling or limiting the spread of disease on this sort of network, and 
so have focussed our attention on edge deletion, which might correspond to forbidden trade pat¬ 
terns or, more reasonably, extra vaccination or disease surveillance along certain trade routes. 
Introducing extra controls of this kind is costly, so it is important to ensure that this is done 
as effectively as possible. Many properties that might be desirable from the point of view of 
restricting the spread of disease can be expressed in terms of forbidden subgraphs: edge-deletion 
to achieve a maximum degree of at most d is equivalent to edge-deletion to a graph avoiding 
the star A'i,d+i, and edge-deletion to maximum component size at most h is equivalent to edge- 
deletion to a graph avoiding all trees on h -|- 1 vertices. Clearly we can also combine criteria of 
this kind, for example edge-deletion to a graph which has maximum component size at most h 
and maximum degree at most d. We are therefore concerned with the following general problem. 

J'-Free Edge Deletion 

Input: A Graph G = (E, E) and an integer k. 

Question: Does there exist E' Q E with \E'\ = k such that G \ E' does not contain any 

F G as an subgraph? 

This is in fact a special case of the more general problem in which we seek to avoid a set 
E of graphs as induced subgraphs (which corresponds to edge-deletion to a hereditary class of 
graphs). We have chosen to focus on the special case of edge-deletion to a monotone class of 
graphs (that is, a class closed under deletion of vertices and edges) as it is reasonable to assume 
in epidemiological applications that if we wish to avoid some subgraph E then we also wish to 
avoid any graph E' obtained from F by adding edges. Moreover, this assumption improves the 
running time of the algorithm by decreasing the size of the family E we wish to avoid, compared 
with expressing our target class in terms of forbidden induced subgraphs (for example, only one 
forbidden subgraph is required to define the class of graphs with maximum degree at most d, but 
to express this in terms of forbidden induced subgraphs we would have to forbid every induced 
subgraph on d-|-2 vertices that contains a vertex of degree d-|-l). However, it is straightforward 
to adapt the algorithm described in Section [3] to consider induced subgraphs. The algorithm 
can also easily be adapted to deal with different costs associated with the deletion of different 
edges (so as to decide whether it is possible to delete edges with a total cost of at most k to 
remove all copies of subgraphs from E). 

A special case of particular interest is the situation in which E is the set of all trees on 
h + 1 vertices, so that we are deleting edges in order to obtain a graph in which every connected 
component contains at most h vertices; h is then an upper bound on the number of vertices 
which may, in the worst case, be infected from a single initially infected vertex. We denote by 
Eh+i the set of all trees on h + 1 vertices, so this special case is the problem 7/i+i-Free Edge 
Deletion. For this case, we also consider two straightforward extensions of the problem that 
are relevant for real-world applications: 

• assigning different weights to different vertices (e.g. corresponding to the number of 
animals in a particular animal holding), and seeking to bound the total weight of each 
connected component; 

• imposing different limits on the size of components containing individual vertices (for 
example, we might want to enforce a smaller size limit for components containing certain 
vertices considered to be particularly high risk). 
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Even the special case T^+i-Free Edge Deletion of our general problem is intractable in 
general for constant h, as demonstrated in the following proposition. The reduction replies on 
the observation that the maximum number of edges in a graph having maximum component 
size h is obtained if the graph is a disjoint union of h-cliques. 

Proposition 1.1. 7h+i-FREE Edge Deletion is NP-complete, even when h = 3. 

Proof. First, observe that an edge set of size k to be deleted suffices as a certificate for this 
problem, therefore 7h+i-FREE Edge Deletion is in NP. We prove NP-hardness by means of 
a reduction from the following problem, shown to be NP-hard in m- 

Perfegt Triangle Cover 

Input: A graph G = {V,E). 

Question: Does there exist a set of vertex-disjoint triangles that cover all vertices in the 

graph? 

Starting with an instance G = (E, E) of Perfegt Triangle Cover (where G has n 
vertices), let G = (E, E) and k = \E\ — n be an instance of TI-Free Edge Deletion; we claim 
that G is a yes instance of Perfegt Triangle Cover if and only if [G, k) is a yes instance 
of 74-Free Edge Deletion. 

Suppose first that G = (E, E) is a yes instance of Perfect Triangle Cover, so there 
exists a set of vertex-disjoint triangles that cover all vertices of G; let E' be the set of edges 
induced by those triangles. Note that E' contains exactly n edges, so the edge-set E\E' is of 
size |E| — n = k] moreover, as there are no edges in E' between distinct triangles, the graph 
G\{E\E') contains no connected component on more than three vertices. Thus {G,k) is a yes 
instance of T^-Free Edge Deletion. 

Conversely, suppose that (G, k) is a yes instance of TI-Free Edge Deletion. Then there 
exists an edge set F such that \F\ = k = \E\ — n and every connected component of G\F is of 
size three or less; by a pigeonhole argument, these components must all be triangles. Therefore 
G is a yes instance of Perfect Triangle Cover, as claimed. □ 

In order to develop useful algorithms for real-world applications, we therefore need to exploit 
structural properties of the input network. In Section [2] we provide evidence that many animal 
trade networks of interest are likely to have small treewidth, and in Section [3] we describe 
an algorithm to solve J^-Free Edge Deletion whose running time on an n-vertex graph of 
treewidth w is bounded by if every graph in E has at most r vertices; this algorithm 

is easily adapted to output an optimal solution. In Section S] we then improve on this for the 
special case of 7h+i-FREE Edge Deletion, describing an algorithm to solve this problem in 
time 0{{wh)‘^^n) on an re-vertex graph of treewidth w. Many problems that are thought to be 
intractable in general are known to admit polynomial-time algorithms when restricted to graphs 
of bounded treewidth, often by means of a dynamic programming strategy similar to that used 
to attack the problem considered here; however, to the best of the authors’ knowledge, the 
usefulness of such algorithms for solving real-world network problems has yet to be investigated 
thoroughly. 

In reality, policy decisions about where to introduce controls are likely to be influenced by a 
range of factors, which cannot all be captured adequately in a network model. Thus, the main 
application of our algorithm will be in comparing any proposed strategy with the theoretical 
optimum: a policy-maker can determine whether there is a solution with the same total cost 
that results in a smaller maximum component size. We provide an example of an experimental 
application of our algorithm to cattle trading networks in Section 14.51 
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In the remainder of this section, we begin by reviewing previous related work in Section 11.11 
before introducing some important notation in Section 11.21 and reviewing the key features of 
tree decompositions in Section [131 

1.1 Review of previous work 

From a combinatorial perspective, we are concerned here with edge-deletion problems. An edge- 
deletion problem asks if there is a set of at most k edges that can be deleted from an input 
graph to produce a graph in some target class. In contrast to the related well-characterised 
vertex-deletion problems [20], there is not yet a complete characterisation of the hardness of 
edge-deletion problems by target graph class. 

Yannakakis jSS] gave early results in edge-deletion problems, showing that edge-deletion 
to planar graphs, outer-planar graphs, line graphs, and transitive digraphs is NP-complete. 
Subsequently, Watanabe, Ae and Nakamura [28] showed that edge-deletion problems are NP- 
complete if the target graph class can be hnitely characterised by 3-connected graphs. There 
are a number of further hardness results known for edge-deletion to well-studied graph classes, 
including for interval and unit interval graphs cographs jH], and threshold graphs [25| and, 
as noted in [ 23 |, hardness of edge-deletion to bipartite graphs follows from the hardness of a 
MAX-CUT problem. Natanzon, Shamir and Sharan |21| further showed NP-completeness of 
edge-deletion to disjoint unions of cliques, and perfect, chain, chordal, split, and asteroidal- 
triple-tree graphs, but also give polynomial-time algorithms, in the special case of the input 
graph having bounded degree, for edge-deletion to chain, split, and threshold graphs. 

Given the large number of hardness results in the literature, it is natural to consider the 
parameterised complexity of these problems. Cai [5] initiated this investigation, showing that 
edge-deletion to a graph class characterisable by a finite set of forbidden induced subgraphs is 
fixed-parameter tractable when parameterised by k (the number of edges to delete): he gave an 
algorithm to solve the problem in time 0(r^^ • where n is the number of vertices in the 

input graph and r is the maximum number of vertices in a forbidden induced subgraph. Further 
fpt-algorithms have been obtained for edge-deletion to split graphs m and to chain, split, 
threshold, and co-trivially perfect graphs m- When considering graphs of small treewidth, 
our algorithm for J-'-Free Edge Deletion (and indeed its adaptation to deal with forbidden 
induced subgraphs) represents a signihcant improvement on Cai’s algorithm, with our running 
time of n. While the fixed parameter tractability of this problem (parameterised by r, 

the maximum number of vertices in any element of J-) restricted to graphs of bounded treewidth 
does follow from the optimization version of Courcelle’s Theorem 13 El, this does not lead to a 
practical algorithm for addressing real-world problems. Note that Proposition 11.11 implies that 
parameterisation by r alone will not be sufficient to give an fpt-algorithm. 

The specific problem of modifying a graph to bound the maximum component size has 
previously been studied both in the setting of epidemiology m and in the study of network 
vulnerability naiHi. The edge-modihcation version we consider here appears in the literature 
under various names, including the component order edge connectivity problem US] and the 
minimum worst contamination problem |21j . Li and Tang |21j show that it is NP-hard to 
approximate the minimisation version of the problem to within 2 — e, while Gross et. al. m 
describe a polynomial-time algorithm to solve the problem when the input graph is a tree. 

1.2 Notation and definitions 

Unless otherwise stated, all graphs are simple, undirected, and loopless. For graph G = {V,E), 
V = V{G) is the vertex set of G, and E = E{G) the edge set of G. We denote the sizes of 
the edge and vertex sets of G as e{G) = |E(G)| and v{G) = |U(G)|. For a vertex v G V{G), 
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we say that vertex u G V{G) is a neighbour of v if {u,v) G E{G), and write Ndv) for the set 
of neighbours of v in G. If [/ C V{G), we write G[U] for the subgraph of G induced by the 
vertex-set U. Given a graph G = {V,E) and a vertex v G V{G), we write G\v for the graph 
G\y \ {u}]. Given a set of edges E' C E{G) and a vertex v G V{G), we write E' \ v for the set 
of edges in E' that are not incident with v. For any vertex v G V{G), we write dciv) for the 
degree of v in G; when the graph G is clear from the context we may omit the subscript. For 
further general graph notation, we direct the reader to m- 

Given two graphs Hi and H 2 with v{Hi) < v{H 2 ), an embedding of Hi into H 2 is an 
injective function 9 : V{Hi) V{H 2 ) such that 6{u)6{v) G E{H 2 ) whenever uv G E{Hi). 
Thus the graph G contains the graph T as a subgraph if and only if there is an embedding 
of E into G. We say that the embedding 0 is a strong embedding (or induced embedding) if 
we have 6{u)6{v) G E{H 2 ) if and only if uv G E{Hi) (so 6 preserves non-adjacency, as well as 
adjacency). 

Given a function 9 which maps a subset of X to a subset of Y, we write Dom(0) and lm(0) 
for the domain and image of 9 respectively. Given a subset X' C Dom(0), we write 9\x' for the 
restriction of 9 to X'. 

A partition P of a set A is a collection of disjoint, non-empty sets whose union is X. We call 
each set in the partition a block of the partition, and every partition corresponds to a unique 
equivalence relation on X where x ~ y if and only if x and y belong to the same block of X. 
If V and V are partitions of A, we say that V refines V if every block of V is contained in 
a single block of P. If "P is a partition of A, and y G A, we write V \y for the partition of 
A \ {y} obtained by removing the occurrence of y from V (and, if this results in an empty set 
in the partition, also removing this empty set). 

1.3 Tree decompositions 

In this section we review the concept of a tree decomposition (introduced by Robertson and 
Seymour in [5B]) and introduce some of the key notation we will use throughout the rest of the 
paper. 

Given any tree T, we will assume that it contains some distinguished vertex r(T), which we 
will call the root of T. For any vertex u G R (T) \ r{T), the parent of v is the neighbour of v on 
the unique path from v to r{T)] the set of children of v is the set of all vertices u G V{T) such 
that V is the parent of u. The leaves of T are the vertices of T whose set of children is empty. 
We say that a vertex u is a descendant of the vertex u if u lies somewhere on the unique path 
from u to r(T) (note therefore that every vertex is a descendant of the root). Additionally, for 
any vertex v, we will denote by the subtree induced by v together with the descendants of v. 

We say that {T,T>) is a tree decomposition of G if T is a tree and T> = {T>{t) ; t G V{T)} is 
a collection of non-empty subsets of V{G) (or bags), indexed by the nodes of T, satisfying: 

1- ^(G) = Ugy(r)^W> 

2. for every e = uv £ E{G), there exists t G V{T) such that u,v £ E>{t), 

3. for every v £ V{G), if T{v) is defined to be the subgraph of T induced by nodes t with 
V £ T>{t), then T{v) is connected. 

The width of the tree decomposition {T,T>) is defined to be \'E{t)\ — 1, and the 

treewidth of G is the minimum width over all tree decompositions of G. 

We will denote by V) the set of vertices in G that occur in bags indexed by the descendants 
of t in T. Thus, V) = [jt'&v{Tt) ^(^0- 

Later in this paper, we will exploit two useful properties that follow from the definition of 
a tree decomposition: 
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1. If u € T){t) and t' is a child of t with v ^ then any path in G from n to a vertex 

w G Vt' must include at least one vertex of 'D{t) \ {n}. 

2. If t G r and ti, t 2 are children of t, then any path in G from a vertex ui G to a vertex 
V 2 G ht 2 must contain at least one vertex from I’(t). 

Although it is NP-hard to determine the treewidth of an arbitrary graph [2] , it is shown in [3] 
that the problem of determining whether a graph has treewidth at most re, and if so computing 
a tree-decomposition of width at most rc, can be solved in linear time for any constant w. 

Theorem 1.2 (Bodlaender [3]). For each w (z N , there exists a linear time algorithm, that 
tests whether a given graph G = {V, E) has treewidth at most w, and if so, outputs a tree 
decomposition of G with treewidth at most w. 

A special kind of tree decomposition, known as a nice tree decomposition, was introduced 
by Kloks m- The nodes in such a decomposition can be partitioned into four types (examples 
in Figure [T]): 

Leaf nodes: t is a leaf in T. 

Introduce nodes: t has one child t', such that Eft') C E{t) and |3^(t)| = \E{t') \ -|- 1. 

Forget nodes: t has one child t', such that E{t') D E{t) and |T’(t)| = \E{t') \ — 1. 

Join nodes: t has two children, ti and t 2 , with E{ti) = E{t 2 ) = Eft). 



Figure 1: The four types of node in a nice tree decomposition. From left to right: a leaf, an 
introduce node, a forget node, and a join node. 

Any tree decomposition can be transformed into a nice tree decomposition in linear time: 

Lemma 1.3 ([E]). For constant k, given a tree decomposition of a graph G of width w and 
0(n) nodes, where n is the number of vertices of G, one can find a nice tree decomposition of 
G of width w and with at most 4n nodes in 0(n) time. 


2 Treewidth of real networks 

While the overall graph of cattle trades in Great Britain from 2001 to 2014 is fairly dense, many 
of the edges are repeated or parallel trades: that is, a farm sending animals over time to the same 
place, or many individual animals being moved at the same time; if we restrict our attention 
to a limited time frame, and ignore movements that would generate multiple edges (that is, we 
require our graph to be simple), the graph is quite sparse: for example, the aggregated graph 
of cattle trades in Scotland in 2009 has an edge-to-vertex ratio of approximately 1.15 when 
aggregating over January alone, and approximately 1.34 when aggregating over the entire year. 
When considering an epidemic, it is much more relevant only to consider trades occurring within 
some restricted time frame (whose precise duration depends on the disease under consideration). 
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Figure 2; A plot of an upper bound treewidth of the largest component in an aggregated, undi¬ 
rected version of the cattle movement graph in Scotland in 2009 over a number of different days 
included: all day sets start on January 1, 2009. Treewidths below eight are exact, treewidths 
over eight are upper bounds of the true treewidth. 


We have investigated the treewidth of real networks arising from two kinds of of cattle 
trade data. First of all, for years between 2009 and 2014, we generated a graph from a type of 
persistent trade link recorded by BCMS in Scotland. The largest of these is derived from the 
trades in 2013, and includes approximately 7,000 nodes and 6,000 edges (this lower density is 
typical when considering only persistent trade links, or trades over a restricted time period). 
None of these graphs has treewidth more than fifteen, with most having treewidth less than 
four. Secondly, in addition to these persistent trade links, we have computed an upper bound 
of the treewidth of the largest component of an aggregated, undirected version of the overall 
network of cattle trades in Scotland in 2009 over a variety of time windows, as illustrated in 
Figure El The treewidths of these components remains low even for large time windows: for 
an aggregation of all movements in a 200-day window the treewidth is below 10, and for all 
movements over the year it is below 18. It is unlikely to be necessary to include a full year of 
movements in the analysis of any single epidemic, as the time scale of most exotic epidemics is 
much shorter. 

When modelling disease processes on graphs, it is common to model contact networks with 
any of a variety of random graph models, including Erdos-Renyi random graphs or random 
graphs with scale-free degree distributions. The cattle trading graphs considered here exhibit 
much lower treewidth than we would expect from an Erdos-Renyi graph with the same density: 
Gao [ID] has shown that Erdos-Renyi graphs with an edge-to-vertex ratio of at least 1.073 
(which our graphs satisfy), have treewidth linear in the number of vertices with high probability. 
Scale-free random graphs have also been shown to have high treewidth: Gao also shows that 
scale-free graphs produced by preferential attachment with at least 11 edges added at each step 
have treewidth linear in the number of vertices with high probability; however, our graphs have 
lower density than this. 

We have excluded markets from the data used above (so that if animals are sold from farm 
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A to farm B via a market, this is considered as a direct trade from A to B), so the low treewidth 
cannot be explained simply by the fact that most trades take place through a relatively small 
number of markets. One possible explanation for the low treewidth is the structure of the 
industry: farms can sometimes be characterised by “type”, with breeders producing calves who 
then might be grown at one or two other farms before eventual slaughter, meaning that cycles 
are unlikely to occur frequently in the network. A more thorough investigation of underlying 
processes that lead to trade networks of low treewidth, and their influence on other graph 
parameters, seems a fruitful direction for future research. 

While we have by no means completed an exhaustive study of the structural properties of 
real-world livestock trade networks, the evidence given here seems sufficient to suggest that 
algorithms which achieve a good running time on graphs of bounded treewidth will be useful 
for this application in practice. Indeed, the benefits of exploiting the treewidth of the input 
graph are demonstrated by our initial experimental results in Section [4.51 

3 The general algorithm 

In this section, we describe an algorithm which, given a graph G together with a nice tree 
decomposition (T, V) of G of width at most re, solves J'-Free Edge Deletion on input G in 
time n. Since there exist linear-time algorithms both to compute a tree-decomposition 

of any graph G of fixed treewidth w, and to transform an arbitrary tree decomposition into a 
nice tree decomposition, this in fact gives an algorithm which takes as input just a graph G of 
treewidth at most w. Thus, we prove the following theorem. 

Theorem 3.1. There exists an algorithm to solve J'-Free Edge Deletion in time on 

an input graph with n vertices whose treewidth is at most w, if every element of J- has at most 
r vertices. 

As with many algorithms that use tree decompositions, our algorithm works by recursively 
carrying out computations for each node of the tree, using the results of the same computation 
carried out on any children of the node in question. In this case, we recursively compute the 
signature of each node: we define the signature of a node in Section f3.11 It is then possible to 
determine whether we have a yes- or no-instance to the problem by examining the signature of 
the root of T. 

We may assume, without loss of generality, that every element of J- contains at least one 
edge: if some F ^ B has no edges, then any graph on at least v{F) vertices will contain a copy 
of F (no matter how many edges we delete) so there cannot be a yes-instance with r or more 
vertices, and a brute-force approach will achieve the desired time bound on input graphs having 
fewer than r vertices. 

In Section 13.21 we describe how we compute the signature of a bag indexed by a given 
node from the signatures of its children, before discussing the running time and a number of 
extensions in Section 13.31 

3.1 The signature of a node 

In this section we describe the information we must compute for each node, and define the signa¬ 
ture of a node. Throughout the algorithm, we need to record the possible states corresponding 
to a given bag. A valid state of a bag T>{t) is a pair consisting of: 


1. a spanning subgraph H of G['D{t)], and 


2. a function (p : Xjr jj —>• {0,1} (where Xjr jj consists of all pairs (F', 6) such that F' is an 
induced subgraph of some F ^ F and 9 is an embedding into H of some induced subgraph 
F" of F') which satisfies the following conditions: 

(a) (p{F',9) = 1 whenever 6 is an embedding of F' into H; 

(b) if (piFi, 9) = 1 and F 2 is an induced subgraph of Fi then 4>{F2,9\v(F2)) = 1; 

(c) if 9) = 1 (where Fi is an induced subgraph oi F £ F and 9 is an embedding of 
F[ into F{) then, if there exist v G V{H) \ lm(0) and u G V{F) \ V{Fi) such that, 
for every w G V{Fi) with uw G E{F), 

• w £ F[, and 

• 9{w)v £ E{H), 

then (t){E\y {Fi) U {u}], 9) = 1, where 9 extends 9 by mapping u to v. 

(d) for every F £ F and every 9, 4>{F, 9) = 0. 


Intuitively, the state of a node tells us which forbidden subgraphs appear in G\yt] (and 

where these appear), so we know which partial forbidden subgraphs we must avoid extending 

to full copies of these graphs. Conditions 2(b) and 2(c) are to ensure consistency, so that if a 

particular partial embedding of a forbidden subgraph is present then the appropriate extensions 

and restrictions of this embedding are also present. 

For any bag P(t), we denote by st{t) the set of possible states of T>{t). Note that there 

are at most 2"^ possible spanning subgraphs H, since \'F>{t)\ < re + 1 (as the treewidth of 
/ 2 

G is w) and 2v 2 < 2^ . For any given spanning subgraph H and any F which is an 

induced subgraph of some F £ F (with v{F") < v{E{)), the number of embeddings of F" 
into F[ is at most the number of injective functions from V{F") to V{E{), which is equal to 
{vFFv{F"))\ ^ < (tc + 1)’'*'^”^; summing over all possibilities for F", we have that the 

total number of pairs {F', 9) such that F” is an induced subgraph of F' and 0 is a an embedding 
of F" into H is at most 


E E 

F ' an induced 
subgraph of F 


v { F ) i 

Y. (»+ir<^"> = EEE 

F^F i=l 3=1 
v { F ) 


F " an induced 
subgraph of F ' 


v{F)\ /i 


( (re + 1)-^ 


v{F) 


SEE 

FgJf i=l 

<\F\{w+^r. 
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Not such pairs will belong to Xf,h, as they may not obey the condition on edges, but this 
certainly provides an upper bound on \Xf,h\- The number of possibilities for the function (p is 
therefore at most Overall, this gives an upper bound of on the number 

of valid states for each bag. Note that, since every element of F has at most r vertices, the 
cardinality of F is at most ^ bound the number of valid states by 

(w+sy 

For any given state a = {H, 9) £ st(t), we set £{t, a) to be the set of edge-sets E' C E{G[Vt]) 
such that Gt = G[Vt] \ E' has the following properties: 

1. Gt[V{t)] = H, 
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2. Gt does not contain any F ^ T as a subgraph, and 

3. for every F ^ F and every subgraph F' of F, if 9 is an embedding of F' into Gt and 
e{V{F'))nV{t) + 0 then c^{F\9\j,^t)) = 1- 

Note that, whenever u is a valid state for t, the set £{t,a) will be non-empty: setting 
E' = E{G[Vt]) will always satisfy all three conditions, since every F & F contains at least one 
edge. 

Since we are interested in determining whether it is possible to delete at most k edges to 
obtain a graph that does not contain any element of as a subgraph, we will primarily be 
interested in a subset of £{t, a): for any node t and a € st(t) we define this subset as 

£k{t, a) = {E' S £{t, a) : \E'\ < k}. 

We then dehne 

dek(t,cj)= min \E'\, 

adopting again the convention that the minimum, taken over an empty set, is equal to infinity. 
Finally, we define the signature of a node t to be the function sig^ : st{t) —>■ {0,1,..., A;, oo} such 
that sig^(cr) = delfc(t,cr). 

Our input graph is then a yes-instance to J^-Free Edge Deletion if and only if there 
exists some a € st(tr) such that sigj^(cj) < k, where tr is the root of the tree indexing the 
decomposition. 

3.2 Computing signatures recursively 

In this section we describe how the signature of any given node in a nice tree decomposition can 
be determined from the signatures of its children (if any). This relationship depends on the type 
of node under consideration, and we will again discuss each type of node in turn, describing 
how to compute delfc(t, a) for an arbitrary state cj. 


Leaf nodes 

Note that, in this case, we have Vt = F>{t), and hence for any state a = a graph Gt 

satisfies the conditions in the definition of £{t, a) if and only H = Gt- Thus 


delfc(t,cj) 


e{Gt) - e{H) 


Introduce nodes 

Suppose that the introduce node t has child t', and that T>{t) \ T>{t') = {u}. Given a state a G 
st(t), we define the introduce-inherited state of t' with respect to a to be the state a' = (i/', (/>'), 
where H' = H \ v and for all {E',9) G Xjr jji we have (l)'{F',6) = (l){E',9). It is clear that a' 
does indeed belong to st(t'). 

We make the following claim concerning the signatures of t and T; recall that dciv) denotes 
the degree of the vertex v in the graph G. 

Lemma 3.2. Let t be an introduce node with child t', and let a = {H,V,c) G st(t). Then 
delk{t,a) = [delfc(t',cr') -h (dcpW] (w) “ > 

where a' is the introduce-inherited state oft' with respect to oj. 
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Proof. Observe first that the edges incident with v in any set E' G £(t, a) must be precisely the 
set Ej; = {uv G G\Vt] : uv ^ E{H)}. We shall in fact demonstrate that 

£{t,a) = {E^UE' : E' G £{t',a')} , 

and hence that 


a) = {^Ey iJ E' ■. E' ^ £k{t', <y'), and \Ey U E'\ < /c| . 

The result will then follow immediately by taking the minimum cardinality of an element in 
this set, using the observation that Ey is disjoint from every E' G £k{t',a'). 

We hrst show that £{t,a) C {Ey U E' : E' £{t',a')}. Let E G £{t,a): we claim that there 
exists a' G iT[t^]intr such that E\Ey G £{t', a'). Set Gt = G[Vt]\E. It suffices to check that Gt[Lt'] 
satisfies the three conditions in the definition of £{t', a'). The first condition is trivially satished, 
by definition of and the second follows immediately from the fact that E G £{t,a) and 

hence Gt does not contain any T G T” as a subgraph. For the third condition, suppose that F' 
is an induced subgraph of some F ^ E and that 9 is an embedding of F' into Gi[I4']- In this 
case 0 is also an embedding of F' into Gt and so (since E G £’(t, a)) we must have 4>{F', 6) = 1; 
hence, by definition of (T[F]intr, we see that (j)'{F',9) = 1, as required. Thus Gt[Lt/] satisfies the 
conditions in the definition of £{t', a') and we conclude that £{t, a) C {Ey U E' ■. E' ^ £{t', o'')}, 
as required. 

Conversely, we now show that {Ey U E' : E' £ £(t',a')} C £{t,(j). Let E' G £{t',a') for 
some a' = G cr[F]intr; we claim that Ey U E' £ £{t,a). Set Gt = G[Vj'] \ E', and let 

Gt = G[Vt]\{E'DEy). It suffices to demonstrate that Gt satishes the conditions in the definition 
of £{t, a). 

The first condition holds by definition. For the second condition, suppose for a contradiction 
that Gt contains a copy of some F £ E; let -0 be an embedding of F into Gt. Since E' £ £{t', a'), 
we know that F is not a subgraph of Gt , so i/j must map some vertex of T to u. It follows that 
'tplv, is an embedding of some F' C F into Gt and hence by dehnition of £{t', a') we must have 
(p'{F', iflv^i) = 1. By definition of (T[t']intr we then have (p{F', iply^) = 1 and hence, by condition 
2(c) in the definition of a valid state, we would also have (p{F,fj) = 1 (since the properties of 
the tree decomposition ensure that every neighbour of 'il)~^{v) in F must belong to £>{£)). It 
then follows that a is not a valid state for t (by condition 2(d) of the definition of a valid state), 
giving the required contradiction. 

For the third condition, suppose that F £ E, F' is an induced subgraph of F, and 9 
is an embedding of F' into Gt. If u ^ ^(-^0 then 9 is an embedding of F' into Gt so, since 
E' £ £{t', a'), we know that 4>'{F', 9) = 1; by the definition of (T[t']intr it follows that 4>{F', 9) = 1, 
as required. Thus we may assume that v £ 9{F')-, suppose that u = 9~^{v). In this case we see 
that 6*|y(F')\{n} is an embedding of F' \u into Gt , so by dehnition of £{t',a') we must have 
(p'{F' \ u, 0|y(F')\{u}) = 1 and hence by dehnition of (T[t']intr we have (p{F' \ u, 0|i7’/\u) = 1. Since 
we know u is a valid state for t, it then follows from condition 2(c) in the dehnition of a valid 
state that we must also have 4>{F',9) = 1, as required. 

This completes the argument that {Ey D E' : E' £ £{t',a'),Ey C Ey} C £(t,a), and hence 
the proof. □ 

Forget nodes 

Suppose that the forget node t has child t', and that £>{£) \ T>{t) = {u}. In this case, we need 
to consider the set of forget-inherited states for t'. Given cr = {H,(j)) £ st{t), we write cr[T]fgt 
for the set of states {H',(p') £ st{t') satisfying: 
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1. H'[V{t)] = H, and 

2. rni) = {(i^',0|p(t)):(i^',0)G W(l)}- 

We make the following claim concerning the signatures of t and t'. 

Lemma 3.3. Let t be a forget node with child t', let a = G st(t), and let (T[t']fgt be the 

set of inherited states oft'. Then 

delfc(f, cj)= min delfc(t', cr'). 

cr'G(T[4']fgt 

Proof. We shall in fact demonstrate that 

£it,cr)= U 

CT'e(T[t']fgt 

and hence that 

£k{t,(T)= IJ {Sk{t',a')}, 

o-'€cr[t']fgt 

from which the result will follow immediately by taking the minimum cardinality of an element 
in each set. 

We begin by showing that Uo-'eo-fi'jfgt ^ S{t,cr). Suppose that E' G £(t',a') for 

some a' G cr[t']fgt; we must verify that E' also satisfies the three conditions in the definition of 
£{t,a). Note that Vt = Vf, and set Gt = G[kt] \ E'. The first two conditions follow trivially 
from the definition of cr[t']fgt and the fact that E' G £{t',a'). For the third condition, observe 
that if 6 is an embedding of E' into Gt then we must have (j)'{E',6\x)(t')) = 1 hence, by 
definition of (T[t']fgt we also have 0(F', 0|x)(i)) = 1, as required. 

Conversely, we now argue that £{t,a) C Let E G £{t,a). We claim 

that there exists a' G o'[T]fgt such that E G £{t',a'). Set Gt = G\yt] \ E = \ E, and let 

a' = where H' = Gt[D{t')] and, for every {F',6) G we have 4>'{F',6) = 1 if and 

only if 4>{F', 6\j)i^t)) = 1- L is straightforward to verify that G u'[t']fgt. We must verify 

that Gt satishes the three conditions in the definition of £{t',a'). The first condition is trivially 
satisfied by our choice of a', and the second follows from the fact that E G £{t,a). For the 
third condition, suppose that 6 is an embedding of F' into Gt] it follows that (p{F', 9\j)(^t)) = 
1 and also that 9\-D{t') is an embedding of F' into El' so, by definition of cf', we also have 
0'(T', 0|x)(i')) = 1, as required. Thus we see that E G £{t',a'), completing the argument that 
£(t,a) C U(T'Go-[t']fgt and hence the proof. □ 


Join nodes 

Suppose that the join node t has children ti and ^ 2 - Note that st(ti) = st(t 2 ) = st{t). In this 
case, we need to consider a set of pairs of inherited states for the two children ti and ^ 2 - We 
write st(ti,f 2 ) for the Cartesian product st(fi) x st(f 2 ) and, given any state a = {E[,(j)) G st(f), 
we write cr[ti, f 2 ]join for the set of pairs ol join-inherited states {{Hi,(j)i),{E[2,(j)2)) £ st(fi,f 2 ) 
satisfying: 

1. Hi=H 2 = H, and 
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2. for every {F',9) G cj) ^(0), and all induced subgraphs Fi and F 2 of F' such that V{Fi) U 
V{F 2 ) = V{F) andV{Fi)nV{F 2 ) = Dom(0), either = 0 or <^2(i^2,0|y(F2)) = 

0 . 

We now make the following claim about the signature of t and those of ti and t 2 - 
Lemma 3.4. Let t be a join node with children ti and t 2 , and let a = € st(t). Then 


delfc(t,cr) 


min {delfc(ti, ai) + delfc(t 2 , 0 - 2 ) - (e(G[2:>(t)]) 

(cri,cr2)€ 

<^[^1 ?^2]join 


eiH)) 


<k 


Proof. We will demonstrate that 

= IJ {EiU E 2 : El € S{ti,ai),E2 e £{t2,a2)}, 

{(ri,a2)£<j[ti,t2] join 


and hence that 

£k{t,a)= IJ {EiU E 2 : El e £k{ti,ai),E2 e £k{t2,(T2),\EiU E2\ < k}. 

(cri,fT2)efT[tl,t2] join 

Notice that, for any Ei € £{ti,ai) and E 2 G £(t2,(72) the properties of the tree decomposition 
ensure that Ei n E 2 = E{G['D{t)]) \ E{H), so the result will follow immediately by taking the 
minimum cardinality of an element from each set. 

We first show that £{t,a) 2 U(ai,t 72 )G<T[ti,t 2 ]join{^i U E 2 : Ei ^ ci ),£"2 G £’(^ 2 , 0 - 2 )}. 
Let ((Ti,(T 2 ) G t 2 ]joiiu and suppose that Ei G £{ti,ai) and £2 £ £’(i2;<72); we claim that 
£1 U £2 G £{t,a). It suffices to show that Gt = G\yt] \ (£1 U £ 2 ) satisfies the three conditions 
in the definition of £{t,a). 

The first condition is trivially satisfied. For the second and third conditions suppose that, 
for some £ G £, £^ is an induced subgraph of £ and 6 is an embedding of F' into Gt- Write 
£/ for F'[Vt^] (for i G {1,2}) and note that 6\vt. is an embedding of £/ into thus, as 

Ei G £{ti,ai) we must have cj)i{F',6\vt.) = 1. Since V{E') = !/(£{) uy(£}), and 9~^{9{V[E'))Pi 
V{t)) = I^(£i) n V{F 2 ), it follows from the definition of <T[ti, t 2 ]join that (/>(£', 0|x)(t)) = 1. This 
demonstrates that the third condition is satisfied. To see that the second condition is also 
satisfied, observe that if £' = £ then we would by the reasoning above have </)(£, ^|xi(t)) = 1, 
contradicting the fact that cj is a valid state, so we see that the second condition must also be 
satisfied. Thus we see that £1 U £2 G £{t, a), as required. 

Conversely, we now show that £{t,a) C U{ai,a 2 )&a[ti,t 2 ]ioiA^^ U £2 : £1 G T(ti,fTi ),£2 G 

£’(^ 2 , < 72 )}- Let E £ £{t,a), and set Gt = G[Vt] \£; we also let £* = £n£(G[Vti]) for i G {1, 2}. 
Now set fjj = {F[,4>i) where 4>i{F',9) = 1 if and only if 9 can be extended to an embedding of 
£' into Gt[Vt^]. 

To verify that (cri,(T 2 ) G <7[ti,t 2 ]join suppose, for a contradiction, that there exists {E',9) G 
and induced subgraphs £1 and £2 of £' such that I^(£i) U V{E 2 ) = I^(£) and I^(£i) H 
V{E 2 ) = Dom(6'), and (l)i{Ei,9\yi^p.^'^) = 4>2{F2-,9\y(^F2)) = 1- Ly definition of 4>i,4>2 it then 
follows (for i G {1,2}) that 0|y(_Fi) can be extended to an embedding 9i of £* into note 

that 9i and 02 agree on V (£ 1 ) H V (£ 2 ) and hence on all vertices of 9~^{T>{t)). We can therefore 
define 9 to be the injective function which agrees with both 9i and 02 on their respective 
domains. It is clear that 0 defines an embedding of F' into Gt and so (as we are assuming that 
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Gt satisfies the conditions in the definition of £{a,t)) we must have = 1; however, 

by definition, 9\x>{t) = 9, so we have 9) = 1, contradicting our initial assumption. We may 
therefore conclude that (iTi,cj 2 ) G t 2 ]join- 

We now proceed to show that Ei G £{ti,ai) for each i G {1,2}. It suffices to demonstrate 
that each \ Ei = Gt[Vt^] satisfies the conditions in the definition of £{ti,ai). The first 

two conditions are trivially satisfied. For the third condition, suppose that F' is an induced 
subgraph of some F £ F and that 9 is an embedding of F' into Gt[Vt^]] it follows immediately 
from the definition of 4 >i that = 1, as required. 

This completes the argument that : Ei G £{ti,ai),E 2 G £’(^ 2 , 0 - 2 )} ^ 

£(t,a), and hence the proof. □ 

3.3 Running time and extensions 

Observe that we can precompute the list of all valid states for a given node t of the tree 
decomposition, as this depends only on the subgraph G['D(t)], for which there are at most 

/ u; + l \ 2 

2v 2 J < 2^" possibilities. To find the set of valid states for a bag which induces a given subgraph 
in G, we first determine all possibilities for FI: this must be a spanning subgraph of G[P(t)], 
and there are most 2'^ possibilities. Next, we determine the set of pairs {F',9) such 

that F' is an induced subgraph of some F £ F and 9 is an embedding into H of some induced 
subgraph F" of F'; to do this, we consider each of the at most + 3)'" possible pairs (F', 9) 

(as described in Section [3.ID . and verify in time at most 0{w‘^) for each such possibility whether 
9 does indeed define an embedding. Finally, for each possible function cj) : — >• {0,1} 

(of which there are at most we must determine whether conditions 2(a)-(d) in the 

definition of a valid state are satisfied: conditions (a) and (d) can each be checked in constant 
time, whereas (b) requires time at most 0(2'") (to consider every possible induced subgraph of 
some Fi G F, which can have at most r vertices) and (c) requires time at most 0{wr'^) (to 
consider every possible choice of a pair of vertices u £ V{F) \ V (Fi) and v £ V {H) \ lm(0), and 
to verify that extending the mapping preserves all edges in F that are incident with u). Thus 
we can generate a lookup table of the valid states for any given node in time 

O ( 2 ^^ ■ 2^^ (|F|(u; + 3)'’ • ^ 

Moreover, we can also precompute, for each possible state a = {H,(f>), the corresponding 
introduce-inherited state and sets of forget-inherited and join-inherited states. 

For an introduce node, there is a single introduce-inherited state for the child node, and this 
can clearly be computed in time 0{\Xjr jji\) = 0{\F{w + 3)*'). 

For a forget node, we need to consider all valid states a' = {H', cj)') for the child. We can first 
verify that H and H' have the correct relationship (taking time at most 0{w‘^)), then we compute 
the set {(F', 0|x)(t)) : {F',9) £ ((/>')~^(1)} (taking time at most 0(|X_F^ff/|) = 0(|F|(rt;-|-3)'’)) and 
verify (in time 0{\Xjr H\)) that this is equal to (/)“^(1). Thus the total time required to compute 
the set of forget-inherited states with respect to a is 0(2"'^+l“^l("'+^)’'t(;^|F|(tc-|-3)^) = 

For a join node, we need to consider all pairs of valid states (ui, 1 T 2 ) for the children, where 
Uj = {Hi,(pi). We can first verify that Hi = H2 = H in time 0{w‘^). For the condition on (p, 
(pi and (p2, we consider each {F',9) £ Xjt^h, determine all possibilities for Fi and F2 so that 
V (Fi) U V (F 2 ) = V (F) (of which there are at most 3^, as each vertex of F' belongs to either Fi, 
F 2 or both) and then verify in constant time whether the condition is satisfied for this choice 
of {F', 9) and (Fi, F 2 ). Thus the total time required to compute the set of join-inherited states 
with respect to a is -|- 3)^3'’) = 

Hence we can perform all precomputation in time Having precomputed the sets 

of states we then, at each of the 0{n) nodes of the nice tree decomposition, iterate over all 


14 


of the valid states for the node, of which there are at most For each possible 

state of a given node, we will need to consider a collection of at most 2^^"^ +l-^l(“'+3)’') inherited 
states (or pairs of states, in the case of a join node) of the child node(s), and to perform various 
constant-time operations for each such state. At each of the 0(n) nodes we therefore do 
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work, so this phase of the computation requires time and the overall time complexity 

of the algorithm is therefore We can bound the size of T by 2^"^, to obtain an upper 

bound on the running time that depends only on n, w and r. 

It is straightforward to adapt the algorithm to solve the related problem in which we wish 
to delete edges to that the resulting graph contains no induced copy of any F £ if: we simply 
consider strong embeddings instead of embeddings, and make analogous changes to the definition 
of a valid state. This does not change the running time as a function of n, w and |J^|, but for 
monotone properties we would typically need to consider a larger set of forbidden induced 
subgraphs than forbidden subgraphs, resulting in an increased running time. 

For simplicity, we have only described the most basic version of the algorithm. However, 
it is straightforward to extend it to deal with a more realistic situation in which the cost of 
deleting different edges may vary: in this case, given a cost function / : E{G) N, we instead 
define delfc(f, a) to be SeeE' /(^)- Additionally, if we wish to output an optimal 

set of edges to delete, we can simply record, for each node t and each state a G st(f), a set of 
edges E' € Sk{t,a) such that \E'\ = delfe(f,(T) (note that in general there may be many such 
optimal sets); computing such a set from the relevant sets for the node’s children requires only 
basic set operations. An element of £{tr,(7), where tj. is the root of the tree decomposition 
and delk{tr,cr) = min^gg^^j^) is then an optimal solution for the problem. Neither of these 
adaptations changes the asymptotic running time of the algorithm. 


4 A specialised algorithm for 7 I+i-Free Edge Deletion 

In this section, we turn our attention to the special case of 7/i+i-Free Edge Deletion and 
describe a more efficient algorithm for this important special case of the problem. Specifically, 
we prove the following theorem. 

Theorem 4.1. There exists an algorithm to solve 7h+i-FREE Edge Deletion in time 0{{wh)‘^'^n) 
on an input graph with n vertices whose treewidth is at most w. 

The general strategy is very similar to that employed in our general algorithm, but in this 
special case we can determine whether our input graph is a yes- or no-instance without using all 
of the information contained in the signature of the root of the tree indexing the decomposition. 
Instead we can record an abridged component-signature for each node, reducing the overall 
running time of the algorithm. 

This relies on two key observations. Firstly, an optimal solution will never delete an edge 
uv if there remains another path from u to v, so it suffices to record which vertices in a bag are 
permitted to belong to the same component, rather than the exact subgraph induced by the 
bag. Secondly, we do not need to know precisely which partial forbidden subgraphs intersect 
the bag and how, rather just the size of the component that contains each vertex of the bag. 

We define the component-signature of a node in Section 14.11 before describing mathemat¬ 
ically how we compute the component-signature of a bag indexed by a given node from the 
signatures of its children in Section 14.21 and discussing the running time and a number of ex¬ 
tensions in Section 14.41 As this algorithm achieves a faster running time and answers a question 
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of particular interest from the point of view of our epidemiological application, the specialised 
algorithm is likely to be of more practical use than the general algorithm described in Section 
El we express the procedure in pseudocode in Section 14.31 and give some initial experimental 
results on its application to cattle trading networks in Section 14.51 

4.1 The component-signature of a node 

In this section, we describe the information we compute for each node, and define the component- 
signature of a node. 

Throughout the algorithm, we need to record the possible component-states corresponding 
to a given bag. A valid component-state of a bag T>{t) is a pair consisting of: 

1. a partition V of 'D{t) into disjoint, non-empty subsets or blocks of size at most h, and 

2. a function c : V ^ [h] such that, for each X ^ V, |X| < c{X). 

We will write u v to indicate that u and v belong to the same block of V. 

Intuitively, V tells us which vertices are allowed to belong to the same component of the 
graph we obtain after deleting edges and c tells us the maximum number of vertices which are 
permitted in components corresponding to a given block of the partition. 

For any bag P(t), we denote by cst(t) the set of possible component-states of T>{t). Note 
that there are at most Byj partitions of a set of size w (where Byj is the Bell number) 
and at most functions from a set of size at most w to [h]; thus the total number of valid 
component-states for T>{t) is at most B^h^ < (wh)'^ (although not all possible combinations of 
a partition and a function will give rise to a valid component-state). 

For any given component-state cr = {V, c) G cst(t), we set T^(t, cr) to be the set of edge-sets 
E' C E{G\Vt]) such that Gt = G[I4] \ E' has the following properties: 

1. for each connected component G oi Gp- 

(a) \V{G)\ < h, and 

(b) if Ct = V{C) n T>(t) / 0, then Gt is contained in a single block Xc of V, 

2. for each block X in V, the total number of vertices in connected components of Gt that 
intersect X is at most c{X). 

Note that, whenever cr is a valid component-state for t, the set £^{t, a) will be non-empty: 
setting E' = i?(G[V)]) will always satisfy both conditions. Since we are interested in determining 
whether it is possible to delete at most k edges to obtain a graph with maximum component 
size h, we will (as in Section [3|) primarily be interested in a subset of £'^{t,a): for any node t 
and cr G cst(t) we define this subset as 

£k{t,a) = {E' G £(t,a) : \E'\ < k}. 

We then define 

cdek(t,cr)= min \E'\, 

adopting the convention that the minimum, taken over an empty set, is equal to infinity. To 
simplify notation, given any a,b G N, we define [o]<f, to be equal to a if a < 6, and equal 
to oo otherwise. Finally, we define the component-signature of a node t to be the function 
csigj : cst(t) —>■ {0,1,. .. , fc, oo} such that csigj(cr) = cdelfc(t, cr). 

Just as in Section [3l our input graph is a yes-instance to T/i+i-Free Edge Deletion if 
and only if there exists some cr G cst(tr) such that csigj^(cr) < k, where E is the root of the tree 
indexing the decomposition. 
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4.2 Computing component-signatures recursively 

In this section we describe how the component-signature of any given node in a nice tree 
decomposition can be determined from the component-signatures of its children (if any). As 
before, this relationship depends on the type of node under consideration, and we discuss each 
type of node in turn, describing how to compute cdelfc(t, a) for an arbitrary component-state a. 

Leaf nodes 

Note that, in this case, we have V) = and hence for any component-state a = {V,c) a 

graph Gt satisfies the conditions in the definition of if and only if the only edges in Gt 

are between vertices that belong to the same block of V. Thus 

cdelfc(t, a) = [\{uv € E {G[T>{t)]) : u ^-p u}|]<^,. 


Introduce nodes 

Suppose that the introduce node t has child t', and that T>{t)\'D{t') = {u}. Given a component- 
state cr = (P, c) € cst(t), where V = {^i,... ,Xr} and v € Xr, we denote by cr[T]c-intr the set 
of introduce-inherited component-states of t'. A component-state {V',c') belongs to cr[T]c-intr 
if and only if it is a valid component-state for £>{£) which additionally satisfies the following 
conditions: 

1. V' = {Xi, ..., Xr-i, hi,, Ts}, with s > 1 and Yi U • • • U Tg = \ {u}, 

2. c'(Xj) = c{Xi) for 1 < i < r — 1, and = c{Xr) — 1. 

We make the following claim concerning the component-signatures of t and t'. 

Lemma 4.2. Let t be an introduce node with child t', and let a = {H,'P,c) G cst(t). Then 


cdelfc(t,cr) 


min cdelk{t', a') -|- \{uv € G[Vt] : u '/'v 'i’}! 

VTcfi^jc-intr 


<k 


Proof. Observe first that the edges incident with v in any set E' G £^{t, cr) must contain the set 
Ey = {uv G G[Vt] : u '/'■p v}. As in the proof of Lemma [321 h suffices to prove that 


£'^{t, a) = lEyUE' : E’ € |J £^{t', a'), Ey C Ey 

We first show that £‘^{t,a) C {Ey U E' : E' ^ cr'), £1^ C Ey}. Let E G 

£'^{t, a): we claim that there exists a' G cr[t']c-intr such that E\v £ £^{t', a'). 

Set Gt = G\yt] \ E, and let C be the component of Gt containing v. Suppose that V = 
Xi,..., Xy, with V £ Xy, and that C\v has connected components Ci,..., Ci. We now define V' 
to be the partition {Yi,..., Yr+^-i} where T) = Xj for 1 < i < r — 1, and = V(Gj)n'D(t') 

for 1 < j < £■ We also dehne c' '.V' ^ [h] by setting 

c(Xj) ifl<i<r — 1 

\V{Gi-y+i)\ ifr<i<r-\-i-l, 
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and then set a' = {V'^c'). It is straightforward to verify that a' G We now claim that 

E\v ^ <y')- To prove this claim, we set Gt = G[Tt'] \ (-S' \ n); we need to verify that Gt 

satisfies the conditions in the definition of Condition 1(a) is trivially satisfied, so we 

consider the remaining conditions. 

For condition 1(b), suppose that C is a connected component of Gt such that Gt = V{G) n 
£>{£) 7 ^ 0. Since G is contained in some connected component G of Gt, and we know that 
every connected component of Gt is contained in a single block of V, it follows that G must be 
contained in a single block Xt of V. There are now two cases to consider, depending on whether 
i = r. If we have i ^ r (implying that n ^ C, so C = C), it follows from the definition of V' 
that G is contained in the block Y) oiV'. If, on the other hand, we have i = r, it follows from 
the dehnition of V' that there is some j G {1,... ,^} such that = G. Thus the condition 

is satished in either case. ^ 

For condition 2, consider a block Y of V'. If Y is also a block of V, then (as Gt is a 
subgraph of Gt) this condition is trivially true, so suppose that this is not the case. Then, by 
construction of V', Y = V{G') nT’(t') for some (maximal) connected component G' of Gt , and 
we have c'{Y) = |y(C")|, so the third condition is indeed satished. 

This completes the argument that Gt satishes the conditions in the dehnition of 
implying that E\v € a') and hence E G {Ey U E' : E' ^ £^{t', a'),Ey C Ey}, 

as required. 

Conversely, we now show that {Ey Li E' : E' £ ,^^^£^(t',a'), Ey C Ey} C £^(t,a). 

Let Ey D Ey, and let E' G £‘^{t', a') for some a' = {V, c!) G (T[t']c-intr (where V' = {Yi,..., Yp], 
with Yi = Xi for 1 < i < r — 1 and Y), U • • • U Y), = Xy \ {n}). We claim that Ey Li E' £ £^{t, a). 

Set Gt = G[V)/] \ E', and let Gt = G[Vt] \ {E' U Ey). It suffices to demonstrate that Gt 
satishes the conditions in the dehnition of £'^{t,a). 

For the hrst condition, let C* be a (maximal) connected component of Gt- If n ^ G, then 
C is a connected component of Gt so we know that |y(C')| < h, and moreover that if Gt = 
V{G) n X>(t) 7 ^ 0 (note that V{G) PI £>{£) = V{G) D T>{t') as v ^ G) then Gt is contained in a 
single block Xy of V and hence a single block of V (as V rehnes V). So it remains to consider 
the case that v £ G. Let Ci,..., be the (maximal) connected components of C'\n. Note that 
it follows from the properties of a tree decomposition that Gt D £>{£) 7 ^ 0 for each i, as there 
must be an edge from v to at least one vertex in each Gf, suppose that ut G Cj nT>(F) for each i. 
Then ui,... ,ui all belong to the same component of H as v, so as cr is a valid component-state 
for t we must have ui,... ,ui G Xy. We know that each Gt is a connected component of Gt 
and so V{Gi) n T>{t) must be contained in a single block of V; since Ui £ Xy, it must be that 
V(Gi) n T>(t') C Yj for some r < j < p. Thus 

v{G)nv{t)= U (<^*nP(0)uM c IJ YjVj{v} = Xy, 

l<i<i 

so V(G) n P(t) is contained in a single block of V, as required. 

For the second condition, let X be a block of V. If v ^ X, then X is also a block of V', and 
moreover the vertices belonging to components of Gt that intersect X are the same as those 
belonging to components of Gt that intersect X, so we are done in this case. Thus we may 
assume that v £ X, and that X \ {n} = Y). U • • • U Yp. Then the number of vertices belonging to 
components of Gt that intersect X is exactly one more than the number of vertices belonging 
to components of Gt that intersect W U • • • U Y^. Thus, the number of vertices belonging to 
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components of Gt that intersect X is at most 


l + Y,c{Yi) = l + c{X)-l = c{X), 

i=r 


as required. 

This completes the argument that {E^UE' : E' G cr'),E^ C E^} C £^{t, a), 

and hence the proof. □ 

Forget nodes 

Suppose that the forget node t has child t', and that T>{t') \ 'D{t) = {n}. In this case, we need to 
consider the set of forget-inherited component-states for t'. Given a = (P,c) G cst(t), we write 
iT[t']c-fgt for the set of component-states ('P^c') G cst(t') satisfying: 

1. V = V' \ {n}, and 

2. for any block F in F, if F \ {n} / 0, we have c'(F) = c(F \ {n}). 

We make the following claim concerning the component-signatures of t and t'. 

Lemma 4.3. Let t be a forget node with child t', let a = (7^, c) G cst(t), and let o'[T]c-fgt be the 
set of inherited component-states oft'. Then 

cdelfc(t, (t)= min cdelkit',a'). 


Proof. As in the proof of Lemma 13.31 it suffices to prove that 

(T^e(T[t^]c-fgt 


It is straightforward to verify that, whenever E' G £'^{t',a') for some o' G we also have 

E' G £''{t,o), implying that Uo-'Go-it']^ fgt ^ £''{t,o). We now argue that we also have 

Y Ucr'e<T[t']c fgt Let E G £^{t,o). We claim that there exists o' G o-[T]c-fgt 

such that E G £^{t',o'). 

Set Gt = G[Vt] \ E = \ E. Let G be the (maximal) connected component of Gt that 

contains v. If C D T’(t) = 0, we set V to be the partition obtained from V by adding one 
additional block containing only n; if G H Pit) 0 then this set of vertices must be contained 
in a single block X oiV, and we set V' to be the partition obtained from V by adding v to the 
block X. 

In the case G D V{t) 0, there is a unique choice of c' : V' ^ [h] which will satisfy the 
dehnition of (T[T]c-fgt: we set c'(A') = c{X') for X' ^ A, and c'(AU{n}) = c{X). If GnVft) = 0, 
we define for any Y gV' 


c'(F) 


c(F) if F G F 
|F(C)| ifF = W. 


We now set o' = {V',c'). It is straightforward to verify that o' G (T[t']c-fgt. It remains to 
demonstrate that E G £'^{t',o'); to do so we will argue that Gt satishes all of the conditions in 
the definition of £"{£, 0 '). Condition 1(a) is trivially satished. 
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For condition 1(b), let C be a (maximal) connected component of Gt, and suppose that 
= V(C) n 'D{t') ^ 0. If C[ = {u} then this condition is trivially satisfied, so we may assume 
that l^(C') n P(t) ^ 0. It then follows from the fact that E G that V{C) D P(t) is 

contained in a single block of V and hence of V'; if we additionally have v G C^, the construction 
of V' ensures that v will also belong to this same block of V. 

For the second condition, let X be a block oiV'. li v ^ X then the condition is immediately 
satisfied due to the conditions for E to belong to so we may assume that u G X. If 

X\u 7^ 0 then, by construction of V', there is some vertex u G such that u and v lie in the 
same component of Gt (note also that rt G X); hence the total number of vertices in connected 
components that intersect X \ u is the same as the number in components that intersect X, 
and so as X \ {u} is a block of V this number must be at most c(X \ {u}) = c'{X). Finally, 
if X = {u}, the only connected component of Gt that intersects X is the component G that 
contains v, and by definition we have c'(X) = |I/(C')|. So the third condition is satisfied in all 
cases. 

Thus we see that E G a'), completing the argument that a) C lJ<T'Gcr[f'] f t 
and hence the proof. □ 


Join nodes 

Suppose that the join node t has children ti and t 2 - Note that cst(ti) = cst(t 2 ) = cst(t). In this 
case, we need to consider a set of pairs of inherited component-states for the two children ti and 
t 2 - We write cst(ti, t 2 ) for the Cartesian product cst(ti) xcst(t 2 ) and, given any component-state 
a = {V,c) G cst(t), we write t 2 ]c-join for the set of pairs of join-inherited component-states 
((Pi, Cl), ( 7 ^ 2 X 2 )) G cst(ti,t 2 ) satisfying: 

1. Pi = P 2 = P, and 

2. for every block X of P, c(X) = ci(X) -|- C 2 (X) — |X|. 

We now make the following claim about the component-signature of t and those of ti and 
t2- 

Lemma 4.4. Let t be a join node with children ti and t 2 , and let a = (P, c) G cst(t). Then 


cdelfc(t,fj) 


min 

(o-i ,0-2)6 

0-[il ,^2]c-join 


{cdelfc(ti,ai) -k cdelfc(t 2 ,o- 2 ) - \{uw G E{G['D{t)]) : u /-p tc}] 


<k 


Proof. As in the proof of Lemma 13.41 it suffices to prove that 

£%t,a)= U {EiUE2:Ei€£%ti,ai),E2e£^{t2,a2)}. 

(cri ,(T2)6cr[tl ,t2]c-join 

To do this, we will exploit some simple observations about V), and Vjj. Note that Vj = 
U Vt 2 , and that by the properties of tree decompositions we have Fi H Vt 2 = ^(0- Thus, for 
any set of vertices U C F, we have 

|C/| = |C/nFil + |C/nF,|-|[/nP(t)|. (i) 

Similarly, for any set of edges F C P(G[F])) we have 

\F\ = |P n P(G[FJ)| + |P n P(G[F,])| - |P n E{G[v{t)])\. ( 2 ) 
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We first show that £^{t,a) C (J 


((Tl ,0'2 )e(T[tl ,t 2 ]c-join 


{El U £"2 : -El € £^{ti,(Ti),E2 G £^{t2,o'2)}- 


Let E G £'^{t,a), and let Gt = G[Vt] \ E. We set ui = (£,ci) and <72 = (£,£ 2 ) where, 
for X £ V and i G {1,2}, Ci{X) is defined to be the total number of vertices in connected 
components of Gt\yt^] that intersect X. It is then straightforward to verify that (<7 i,(T 2) G 
i 2 ]c-join- Moreover, if we set .Ei = .E D and E 2 = E (1 then it follows easily 

that -El G £^{ti,ai) and E 2 G £^{t 2 ,<y 2 ). Since it is clear that -E = -Ei U E 2 , this shows that 
£ U(ai,,T2)Ga[ti,t2]c-join{^i U -E 2 : -El G £‘^{ti,ai),E2 G E'=(t2, ^^ 2 )}. 

Conversely, we will now demonstrate that o- 2 )Go-[ti ^ -Ei G E'^(ti, fJi),-E 2 G 

E'^(t 2 )< 72 )} C £^{t,a). To do this, fix {(7i,a2) G (7[ti,f 2 ]c-join) and let -Ei G £^{ti,ai) and 
-E 2 G £^(t2,(72)] setting Gt = G[Vt] \(-Ei U-E 2 ), we then need to demonstrate that Gt satisfies the 
three conditions in the definition of £^{t,a). We also set Gi = \-Ei and G 2 = \ E 2 . 

Observe that, if C is a connected component of Gt with C H £(1) = 0 then G is contained 
entirely in either Gi or G 2 , and so \C\ < h. This fact, combined with the fact (demonstrated 
below) that condition 2 is satisfied, shows that condition 1(a) is satisfied. 

For condition 1(b), consider a connected component G of Gt, where Gt = V{C) PI X’(f) 7 ^ 0. 
We will argue that, for any two vertices n, in G V{G), we must have u and w belonging to the 
same block of V. Recall that ~'p defines an equivalence relation on T’(f). Whenever there is a 
path from ui to U 2 in Gi, ui and U 2 belong to the same component of Gi and hence ui 1 x 2 ; 
similarly, if there is a path from ui to 1 x 2 hi G 2 , we have xxi U 2 - Now consider u,w € V{G). 
Since u and w belong to the same connected component of Gt, there is a path P from xx to tc in 
Gt- Let xxi = XX, 1 X 2 ,..., XX,. = xc be the vertices of V(P) n V{t), listed in the order in which they 
occur as P is traversed from xx to w. Recall from the properties of a tree decomposition that 
any path from a vertex in l/)i \ ^(^) to a vertex in Vt 2 \ Pit) must pass through Pit). Thus, for 
l<x<r — l,it follows that the segment of P from xx, to xx,+i is entirely contained in either Gi 
or G 2 , and hence that ut ~'p xx,+i. By transitivity, this implies that xx ~-p w, or in other words 
that XX and w belong to the same block of P, as required. 

For the second condition, let X be a block of P, and let Gi,... ,Gr be the components of 
Gt that intersect X. We want to show that 


Y,\ViC^)\<ciX). 

i=l 


For each i, set g\^^ = V (G,) fl V), and G^-^^ = V (G,) D Vt 2 ■ Notice that {Gi D X,..., G^ n X} is 
a partition of X, and moreover that G, n X = G^^^^ fl Hence 


Y.\ViGt)\ = Y,\cf''^cf\ 

i=l 

= E(|c"l + |cPl-|c«ncf» 


2=1 


2 = 1 


EiG’i + Eicf’i-EiG'nc.' 


(2) I 


■y(I) o /^(^) I 


2 = 1 
r 


2 = 1 


2 = 1 


EiG’i + EgFi-w. 

2 = 1 2=1 


Moreover, it is clear that Ui=i is the set of vertices in components of Gi that intersect 
X, so IU=1 < ci(X); since the sets G^^\...,Gr^^ are clearly disjoint, this implies that 
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< ci(X). Similarly, we have Yll=i ^ C 2 (X). Thus we see that 

r 

< ci{X) + C 2 iX) + |X| < c{X), 

i=l 

by definition of cr[ti, t 2 ]c-join 5 so the third condition is satisfied. 

This completes the argument that U(o-i (T 2 )ecr[ti t 2 ]c U E 2 : Ei £ £‘^{ti,ai), E 2 G 

< 72 )} C and hence the proof. □ 

4.3 Algorithmic formulation 

In this section, we have expressed our lemmas for recursively calculating component-signatures 
at nodes in a tree decomposition as pseudocode (Algorithms [1] to [5]). For convenience, “infinity” 
is a single arbitrarily large number; in an implementation this might be the maximum number 
possible in the programming system. 


Algorithm 1 Algorithm generating the set of possible component-states st{t) for bag V^t) 

Input: A node t of the nice tree decomposition T, the bag at that node 'D{t), the graph G, 
integer h 

Output: A set of component-states st{t) 


states ■<— 0 

allPartition ■£- all partitions of T’(t) such that each block is of size at most h 
for V G allPartition do 

allFunctions all functions c from V = {Xi,... A|-p|} to [h] such that c{X) > |A| 
for c G allFunctions do 
add {V, c) to states 

end for 
end for 
return states 


Algorithm 2 Algorithm for finding component-signatures of leaf nodes 

Input: A leaf node t of the nice tree decomposition T, the bag at that node P{t), the graph 
G 

Output: A mapping from component-states to del values for component-states of the node 


delValues empty dictionary 

allStates ■£- get all the valid {V,c) component-states using Algorithm [T] 
for {V, c) in allStates do 

if llttrc G e{G['D{t)]) : u '/"p u;}| < k then 

delValues[t, {V,c)] <— \{uw G e{G['D{t)]) : u 9 ^-p rt;}| 
else 

delValues[t, {V,c)] £- infinity 

end if 
end for 

return delValues 
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Algorithm 3 Algorithm for finding component-signatures of introduce nodes 

Input: An introduce node t of the nice tree decomposition T, the bag at that node the 
graph G, the child t' of t, and a table delValuesChild of cdel^ values for component-states 
of t' 

Output: A mapping from component-states to cdel^ values for component-states of the node 


delValues empty dictionary 

allStates ^ get all the valid {V,c) component-states using Algorithm [T] 
for {V = {^ 1 ,... ,Xr},c) in allStates do 

]intr ^ 0 

refinements ■<— all partitions of Xr\v 
for {li,..., y^} G refinements do 
V ^{Xi,...,Xr-l,Yi,...,Ys} 

allCs ■(— set of all functions c' : V ^ [h] such that c'{Xi) = c{Xi) ifl<i<r— 1 and 
ELlC'{Yi)=c{Xr)-l 
for d G allCs do 

add {V\c!) to cr[t']intr 

end for 
end for 

minValue ■<— inhnity 
for a' G (j[t']intr do 

value ■(— delValuesChild[t', a'] -|- \{uv G E{G['D{t)]) : u /-p n}| 
if value < minValue then 
minValue ^ value 

end if 
end for 

if minValue < k then 

delValues[t, (P, c)] ^minValue 
else 

delValues[t, {V,c)] ^ infinity 

end if 
end for 

return delValues 
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Algorithm 4 Algorithm for finding component-signatures of forget nodes 

Input: A forget node t of the nice tree decomposition T, the bag at that node T>{t), the 
graph G, the child t' of t, and a table delValuesChild of cdel^ values for component-states 
of t' 

Output: A mapping from component-states to cdeh values for component-states of the node 


delValues ^ empty dictionary 

allStates •(— get all the valid {V,c) component-states using Algorithm [T] 
for {V, c) in allStates do 

allPartitions ^ all partitions V' such that V = 'P'\{v} 
for V G allPartitions do 
allc' 0 

c' empty function 
vSingleton false 

for Y do 

if 7 ^ 0 then 

c'iY) ^ c(n{^}) 

else 

vSingleton true 

end if 
end for 

if not vSingleton then 
add d to alld 
else 

for i = 1 to h do 
c'({u}) ^ i 
add d to alld 

end for 
end if 

for V € allPartitions do 
for d G alld do 

add {V',d) to CF[t']fgt 

end for 
end for 
end for 

minValue infinity 
for a' G o'[t']fgt do 

value ■(— delValuesChild[t', a'] 
if value < minValue then 
minValue ^ value 
end if 
end for 

if minValue < k then 

delValues[t, (P, c)] ^minValue 
else 

delValues[t, {V, c)] infinity 

end if 
end for 

return delValues 
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Algorithm 5 Algorithm for finding component-signatures of join nodes 

Input: A join node t of the nice tree decomposition T, the bag at that node the graph 
G, the join node’s children and a table delValuesChild of cdel^ values for component- 

states of ti and t 2 

Output: A mapping from component-states to cdelfc values for component-states of the node 


delValues ^ empty dictionary 

allStates get all the valid {V,c) component-states using Algorithm [T] 
for {V, c) in allStates do 
^\tlit‘2\join ^ 0 

1-^2 ^ 

allFunctions ■(— set of all function pairs (ci,C 2 ) such that for every block X £V, c{X) = 
ci{X) + C2{X) - \X\ 

for (ci,C 2 ) € allFunctions do 

add {{Vl,Ci),{V 2 ,C 2 )) to Cr[ti,t 2 ]join 

end for 

min Value ^ infinity 

for (fTi,cJ 2 ) G cr[ti,t 2 ]join do 

value ■(— delValuesChild[ti, cJi] -|- delValuesChild[t 2 , o' 2 ] - \{uw G E{G[V{t)]) : u 9 ^-p 

u;}| 

if value < min Value then 
minValue •(— value 

end if 
end for 

if minValue < k then 

delValues[t, {V, c)] •<— minValue 
else 

delValues[t, {V, c)] ^ infinity 

end if 
end for 

return delValues 
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4.4 Running time and extensions 

At each of the 0{n) nodes of the nice tree decomposition, we will generate, and then iter¬ 
ate over, fewer than (wh)'^ component-states for that node. For each of those component- 
states, we will need to consider a collection of at most (wh)^ inherited component-states (or 
pairs of component-states, in the case of a join node) of the child node(s). In the algorithm, 
we first generate each of the component-states for a given node, and the corresponding set 
of inherited component-states for its children, then iterate over each relevant combination of 
component-states, performing various constant-time operations. Thus, at each of 0(n) nodes 
we do O work, giving an overall time complexity of 0{{wh)‘^'^n). 

For simplicity, we have only described the most basic version of the algorithm; however, it 
is straightforward to extend it to deal with more complicated situations, involving any or all of 
the following. 

Deleting edges so that the sum of weights of vertices in any component is at most h, 
where a weight function w : V{G) — )• N is given: change condition 1(a) in the definition 
of £^{t, a) to ^(^) ^ h, and add to the definition of the set of valid component-states 

for a node the condition that, for each block X of V, we have ^ c(Ai). 

Deleting edges so that each vertex v belongs to a component containing at most 
i{v) vertices, where a limit fnnction i : V{G) — )• N is given: change condition 1(a) 
in the definition of to |F(C')| < min„gy(c') £(u), and add to the definition of the set 

of valid component-states for a node the condition that, for each block X of V, we have 
c{X) < min^ 

Neither of these adaptations changes the asymptotic running time of the algorithm. It is also 
straightforward to deal with different deletion costs for different edges, and to output an optimal 
set of edges to delete, as discussed for the general algorithm in Section 13.31 

4.5 Experimental results 

As an example of its practical use, we have tested an implementation of our algorithm and a 
constraint satisfaction programming (CP) formulation of the problem to find minimum deletions 
to maximum components of five vertices on graphs derived from the persistent Scottish trade 
links. 

We report early preliminary results in Table [H giving the minimum deletions found by our 
algorithm and by the CP, and the time required. For the CP, most of the deletions were not 
confirmed to be minimum within 10 hours, but we report the time to this confirmation if it 
occurred within the 10 hour running time limit. We use the MiniZinc [25] modelling language 
as a front-end to the Gecode solver [211 Ej for our CP. 

Neither our algorithm nor the CP solver were successful in finding a minimum deletion 
within 24 hours on several of the graphs with higher treewidths (15 and 12), and we do not 
report these in the table. In all cases our algorithm is much faster in clock-time than the CP 
approach: both could likely be improved by careful optimisation of the implementations. 

The implementation of our algorithm takes as input a tree decomposition (generated using 
LibTV^fl [7]) of a graph and the graph itself. Our implementation of the tree decomposition- 
based method can be found at github.com/magicicada/fpt-edge-deletion. 

^available at treewidth.com/treewidth/index.html 
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Table 1: Experimental results for deletion to a maximum component size of 5 comparing the 
performance of our tree decomposition based algorithm and the CP method. v{G) is the number 
of vertices in the graph, e{G) the number of edges, and tw{G) the treewidth. For our tree 
decomposition-based method we report the minimum deletion found and the approximate clock 
time to find it. For the CP method we report the minimum deletion found within 10 hours, the 
approximate clock time to hnd it, and the approximate clock time to confirm that deletion is 
minimum: a dash in this column indicates that the confirmation that the deletion was minimum 
did not complete within 10 hours. The graphs used are a selection of anonymised cattle trading 
graphs from Scotland. Times are given in seconds, with times over an hour given to the nearest 
half hour. 


Graph information Tree decomposition method CP method 


Graph 

ID 

v(G) 

e(G) 

tw(G) 

Minimum Time 

deletion 

found 

Minimum Time 

deletion 

found 

Time to 
conhrma- 

tion 

2010-0 

104 

no 

4 

38 

53 

90 

27000 

- 

2010-3 

45 

45 

3 

11 

4 

25 

21600 

- 

2010-4 

38 

38 

3 

20 

7 

24 

25200 

- 

2010-5 

37 

40 

4 

7 

64 

7 

1032 

1046 

2012-0 

97 

119 

5 

58 

5400 

93 

7200 

- 

2012-1 

72 

74 

3 

20 

11 

38 

32400 

- 

2012-4 

31 

30 

2 

12 

3 

12 

1082 

- 

2012-5 

49 

52 

3 

15 

7 

30 

19800 

- 

2013-1 

45 

47 

4 

20 

35 

30 

19800 

- 

2013-3 

61 

62 

4 

19 

38 

36 

7200 

- 

2013-4 

35 

38 

4 

15 

5 

20 

3600 

- 

2013-6 

39 

41 

3 

14 

4 

17 

7200 

- 

2014-0 

32 

49 

4 

28 

445 

28 

2738 

- 

2014-1 

47 

47 

3 

18 

6 

24 

7200 

- 

2014-2 

57 

59 

4 

17 

10 

33 

19800 

- 

2014-3 

41 

40 

2 

15 

5 

31 

360 

- 

2014-5 

31 

35 

4 

8 

18 

8 

412 

476 

2014-6 

48 

49 

3 

21 

11 

21 

10800 

- 

2014-7 

31 

32 

3 

16 

4 

16 

7 

- 
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5 Conclusions and open problems 


We have investigated the relevance of the well-studied graph parameter treewidth to the struc¬ 
ture of real-world animal trade networks, and have provided evidence that this parameter is 
likely to be small for many networks of interest for epidemiological applications. Motivated 
by this observation, we have derived an algorithm to solve J^-Free Edge Deletion on input 
graphs having n vertices and treewidth bounded by some fixed constant w in time 
if no graph in has more than h vertices. The special case of this problem in which J-” is the 
set of all trees on at most h + 1 vertices is of particular interest from the point of view of the 
control of disease in livestock, and we have derived an improved algorithm for this special case, 
running in time 0{{wh)‘^^n). It is straightforward to adapt both algorithms to deal with more 
complicated situations likely to arise in the application. 

Many open questions remain concerning the complexity of this problem more generally, as 
we are far from having a complete complexity classification. We know that useful structure in 
the input graph is required to give an fpt-algorithm: we demonstrated that it is not sufficient 
to parameterise by the maximum component size h alone (unless P=NP). However, it remains 
open whether the problem might belong to FPT when parameterised only by the treewidth 
w] we conjecture that treewidth alone is not enough, and that the problem is W[l]-hard with 
respect to this parameterisation. Considering other potentially useful structural properties of 
input graphs, one question of particular relevance to epidemiology would be the complexity of 
the problem on planar graphs: this would be relevant for considering the spread of a disease 
based on the geographic location of animal holdings (in situations where a disease is likely to 
be transmitted between animals in adjacent fields). 

Furthermore, animal movement networks can capture more information on real-world ac¬ 
tivity when considered as directed graphs, and the natural generalisation of T^+i-Free Edge 
Deletion to directed graphs in this context would be to consider whether it is possible to 
delete at most k edges from a given directed graph so that the maximum number of vertices 
reachable from any given starting vertex is at most h. Exploiting information on the direction 
of movements might allow more efficient algorithms for this problem when the underlying undi¬ 
rected graph does not have very low treewidth; a natural Hrst question would be to consider 
whether there exists an efficient algorithm to solve this problem on directed acyclic graphs. 
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